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I.  INTRODUCTION:  SUMMARY  OF  RESULTS 
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In  recent  years,  it  has  become  clear  that  communication  time  between 
processors  is  one  of  the  most  severe  limiting  factors  in  designing  high 
speed  parallel  computers.  This  being  the  case,  it  becomes  necessary  to 
design  networks  of  processors  in  which  the  communication  time  is  small. 

A  mathematical  version  of  this  design  problem  is  the  problem  of  designing 
graphs  of  given  fixed  degree  and  number  of  vertices  with  small  diameter. 
(The  report  contains  definitions  of  these  terms.)  In  this  report,  we  have 
attacked  this  problem  from  two  directions.  First,  we  have  constructed  a 
heuristic  al gori thm  whi ch  finds  graphs  with  small  diameter,  and  implemented 
it  on  a  computer.  The  program  is  written  in  the  C-language.  We  have 
compared  the  results  with  previous  best  known  results.  Secondly,  we  have 
studied  a  collection  of  graphs  which  have  compact  and  systematic  descrip¬ 
tions,  the  so-called  Cayley  graphs.  These  graphs  have  the  advantage  that 
routing  algorithms  for  them  are  easy  to  specify.  We  have  given  a  crude 
comparison  of  their  diameter .with  that  of  a  theoretical  bound,  and  studied 
a  specific  class  of  them,  the  "modified  cube-connected  cycles." 

The  results  we  have  obtained  may  be  summarized  as  follows. 


(a)  The  heuristic  algorithm  we  have  constructed  is  an  improvement  over 
all  previous  algorithms  of  its  kind.  In  particular,  we  have  improved 
many  of  the  best  known  values  for  dense  graphs  of  given  degree  and 
diameter.  A  complete  description  is  given  in  section  II-B,  where 
Fig.  1  shows  all  the  improved  values  we  have  obtained. 


(b)  Our  algorithm  does  not  find  the  densest  known  graphs  in  cases  where 
they  are  constructed  using  systematic  combinatorial  constructions. 
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This  suggests  that  one  should  study  a  restricted  class  of  graphs, 
having  systematic  descriptions. 

(c)  The  symmetric  groups  Sn  admit  graph  structures  whose  diameters 
approach  the  theoretical  bound  (Moore  bound)  arbitrarily  well  as  n  ». 
This  suggests  that  they  should  be  studied  much  more  carefully  as  a 
possible  source  of  efficient  communications  networks* 

(d)  A  slight  modification  of'the  cube-connected  cycles  of  [P-V]  produces 
an  infinite  family  of  graphs  whose  diameter  grows  as  2  log2(K),  where 
K.  is  the  number  of  points  in  the  graph.  This  compares  favorably  with 
5/2  log2(K),  which  is  the  diameter  of  the  cube-connected  cycles. 

(e)  A  layout  is  given  for  these  modified  cube  connected  cycles,  whose  area 

is  3/2  times  the  area  of  the  cube  connected  cycles  with  the  same 

? 

number  of  nodes.  The  VLSI  measure  of  complexity,  AT  ,  is  thus  in 

fact  slightly  decreased,  by  a  factor  of  24/25. 

> 

The  substantiation  for  these  results  is  obtained  by  running  the 
heuristic  program  which  we  have  designed,  and  by  theoretical  analysis 
contained  in  section  II. 
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II.  DETAILED  DESCRIPTION  OF  RESULTS 
(A)  Measures  of  Communication  Time 

Throughout,  we  will  be  interested  in  arrays  of  "processors,"  connected 
by  "wires."  The  nature  of  these  processors  is  not  specified,  because  we 
want  to  study  the  general  problem  of  communication  time,  without  restricting 
to  a  specific  situation.  We  will  formalize  this  notion  by  considering 
graphs  F  ,  where  the  vertices  of  the  graph  correspond  to  processors  and 
edges  to  wires.  W-2  will  suppose  that  the  array  functions  is  such  a  way 
that  at  every  time  information  is  allowed  to  flow  along  one  wire.  The  time 
required  to  move  along  one  wire  is  presumed  to  be  constant.  By  the  distance 
between  two  processors,  or  the  corresponding  vertices  x  and  y,  we  mean 
the  length  of  the  shortest  path  in  F  from  x  to  y.  A  path  of  length 
n  from  x  to  y  is  a  sequence  of  vertices  <j>  =  (vQ ,  . . . ,  vn>  in  r  , 
so  that  vQ  =  x,  vn  =  y,  and  so  that  for  each  i,  v^v^  is  an  edge  of  T  . 
Clearly,  this  is  the  same  as  the  length  of  the  shortest  sequence  of  wires 
connecting  the  processors  corresponding  to  x  and  y.  We  denote  the 
distance  by  d(x,y).  By  the  diameter  of  the  array,  we  mean 

max  d(x,y)  =  A(x,y)  . 

x.yer 

where  the  max  is  taken  over  all  pairs  of  vertices  in  r  .  (From  this 
point  on,  we  no  longer  speak  of  the  arrays  of  processors  but  only  of  their 
corresponding  graphs.)  We  wish  to  study  the  problem  of  designing  graphs 
with  a  given  number  of  vertices,  having  small  diameter.  Of  course,  with 
no  constraints  on  the  graph,  this  is  a  trivial  problem,  since  complete 
graphs  all  have  diameter  1  .  However,  technology  dictates  that  the  number 


of  edges  from  each  vertex  should  be  less  than  or  equal  to  some  finite  number 
d  .  Accordingly,  we  say  that  a  graph  is  regular  of  degree  d  if  every 
vertex  has  precisely  d  edges  coming  from  it,  and  we  attempt  to  solve  the 


4 


problem  of  minimizing  the  diameter  of  regular  graphs  of  degree  d  having, 
say,  N  vertices.  There  is  an  obvious  a  priori  upper  bound  to  N,  given  d, 
called  the  Moore  bound  (see  [Bo]),  which  is 

N  <  1  +  d  +  d(d-l)  +  ...  +  d(d-l)k-1 

where  k  is  the  diameter  of  the  graph.  This  says  that  asymptotically,  k 
grows  at  least  as  fast  as 


logd-l  ^N)  =  Infer  ' 

It  is  known,  however,  that  ^ bound  is  only  sharp  in  a  finite  number  of 


cases 


for  d  >  3  (for  d  =  2,  the  cyclic  graphs  are  all  examples  where  it  is  sharp), 
and  it  seems  to  be  the  case  that  it  in  fact  is  rather  crude.  Our  efforts 
toward  studying  this  problem  consist  of  the  construction  of  a  heuristic 
algorithm  (see  §B),  and  some  specific  constructions  derived  from  group 
theory  (see  §0). 

The  diameter  is  itself  only  a  weak  measure  of  the  effectiveness  of  the 
processor  array.  Suppose,  for  instance,  that  the  processors  have  no  memory, 
and  that  one  wishes  to  simultaneously  transfer  information  from  the  vertex 
v.  to  the  vertex  w^ ,  for  i  =  1,  ...»  k.  Since  there  is  no  memory, 
one  must  produce  a  collection  of  paths  4>^  of  length  n  in  r  ,  with 
v(i)  y(j)  .for  all  l  <  i,  j  <  k,  and  so  that  vi 1  ^  vj;^  = 


w. 

i 
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Here,  <J>^  =  jvg1^ . vn^  )  *  Such  a  co^ectlon  °f  paths  is  called  a 

k-separated  multipath  in  r  from  (v-| ,  ....  v^)  to  (w^,  ...»  w^).  We 

define  the  k-separated  distance  between  Cv-| ,  ....  v^)  and  (w^ . w^) 

to  be  the  minimum  length  of  all  k-separated  multipaths  from  Cv^ ,  . ..,  vk) 

to  Cw-|,  ...»  w^)  in  T  ,  and  denote  it  hy  d(<(v1  ,  ....  v^) ,  (.w-j . w^). 

The  k-separated  diameter  is 

max  dk  |  ( V'j ,  (w1 ,  ....  wk)J  =  Ak(r) 

where  the  max  is  taken  over  all  pairs  of  k- triples,  with  v.  t  v.,  w.  f  w. 

*  J  *  J 

unless  i  =  j. 

The  k-separated  diameter  Ak(T)  takes  congestion  into  account,  and  so 
is  a  more  sensitive  measure  of  effectiveness  of  the  processor  array. 

However,  it  assumes  that  the  processors  have  no  memory.  We  will  define 
another  measure  of  efficiency,  Ak(F,£),  which  assumes  that  each  processor 
has  £  units  of  memory.  By  a  (k,£)-separated  multipath  of  length  n  from 
(v-j ,  ... ,  vk)  to  (w  -j,  ....  wk) ,  we  mean  a  k-tuple  <j>^  of  paths  of 
length  n  ,  so  that 

(a)  is  a  path  from  v^  to  w. . 

and 

(b)  Each  of  the  k- tuples  IvP^,  ...»  I  contains  each 

vertex  at  most  £  times. 

Note  that  Ak(r,l)  =  Ak(r),  and  that  for  £  >  k,  Ak(r,£)  =  A(r). 

Intuitively,  Ak(r,£)  measures  the  time  required  to  transfer  the  information 
in  any  k-tuple  of  processors  to  any  other  k-triple  of  processors,  given  that 
each  processor  has  £  units  of  memory. 
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We  will  now  observe  that  for  any  graph  r,  A^r.Jl)  is  just  the 

diameter  of  a  graph  associated  with  r  .  For,  form  the  k-fold  product  graph 

Tx  ...  xr  .  In  rx  ...  xr  ,  consider  the  full  subgraph  on  the  set 
0 

£  Vx  ...  xV,  (here  V  denotes  the  vertex  set  of  r),  where 

Vk  =  !<V  *'•’  V 

0 

easy  to  see  that  Ak(r,£)  =  A( ) .  Thus,  we  have  reduced  these  more 
sophisticated  measures  of  effectiveness  to  a  diameter  question;  this  will 
be  a  useful  reduction  in  view  of  the  algorithm  to  be  considered  in  the  next 
section. 

The  diameter  is  unfortunately  often  rather  expensive  to  compute. 
Consequently,  one  would  like  to  obtain  some  less  sensitive,  more  easily 
computable  invariants  of  graphs  which  still  have  some  relation  to  the 
diameter. 

Definition  The  qirth  of  a  graph  r  is  the  length  of  the  shortest  cycle  of  r. 


no  vertex  appears  more  than  l  times  j  .  Then  it  is 


We  note  that  if  a  graph  has  diameter  k,  then  its  girth  is  at  most 
2k+l .  Intuitively,  the  girth  is  inversely  related  to  the  diameter.  For  a 
fixed  number  of  points,  small  diameter  tends  to  imply  large  girth.  We 
summarize  the  facts  known  about  girth  relating  to  this  problem. 


(a)  There  is  a  lower  bound  for  the  number  of  points  in  a  graph  with  a 

given  girth,  analogous  to  the  Moore  bound  (see  [B]). 

(b)  For  a  graph  of'eTd  girth  g  =  2k  +  1 ,  and  degree  d,  the  number  of 

points  is  at  least 


1  +  d  +  d(d-l)  +  ...  +  d(d-l ) 


This  bound  is  attained  only  for  d  =  2,  or  for  d  =  3,  7,  and 
possibly  5,  7,  g  =  5. 
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(c)  For  a  graph  of  even  girth  g  =  2k  and  degree  d,  the  lower  bound  is 

1  +  d  +  d(d-l)  +  ...  +  d(d-l)k“2  +  (d-l)k_1 

This  bound  is  known  to  be  attainable  only  for  g  =  4,  6,  8,  or  12. 

Also,  only  d  =  ps  +  1  are  known  to  occur,  where  p  is  a  prime.  The 
diameter  is  in  each  case  k.  These  graphs  provide  by  far  the  optimal 
graphs  for  the  diameter  problem  with  given  diameter  and  degree. 

We  define  certain  numbers  related  to  the  girth.  Given  a  vertex  x  (  r, 

we  define  Nk(x)  to  be  the  number  of  vertices  connected  to  x  by  a  path 

of  length  k,  and  define  v.(r)  to  be  min  N.  (x).  Note  that  if  T  is  a 

x  €T 

regular  graph  of  degree  d,  then  vk(r)  is  bounded  above  by 
1  +  d  +  d(d-l)  +  ...  +  d(d-l)k_1  =  UjJd),  and  that  the  girth  of  T  is  >  % 
if  vk(r)  =  yk(d)  for  all  k  <  1/2.  Consequently,  to  maximize  girth,  one 
should  attempt  to  maximize  successively  v2(r),  v3(r),  ....  ( r ) ,  vk+1(r), 

...  .  In  each  case,  vk+^(r)  should  be  maximized  subject  to  the  constraint 
that  one  remains  at  an  optimum  for  the  previous  values  of  the  subscript. 
vk  is  quickly  computed  for  small  values  of  k,  and  the  v£s  are  another 
more  computable  collection  of  invariants  of  graphs,  which  are  related  to 
the  efficiency  of  the  associated  array  of  processors.  This  is  particularly 
useful  in  attempting  to  work  with- the  measures  Ak(r,£),  since  the  graphs 

l 

rk  are  usually  quite  large,  so  the  time  spent  computing  invariants  is  of 
primary  importance. 
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(B)  A  Heuristic  Algorithm 

In  this  section,  we  produce  a  "hill-climbing"  algorithm  using  a 
particular  heuristic  criterion  to  find  graphs  of  a  given  fixed  degree  and 
find/number  of  points  with  small  diameter. 

Let  r  be  a  graph  of  degree  d,  and  let  vi »  v2»  wi »  w2  vertlces 
of  T  so  that  v-|  V2  and  w^w 2  are  edges  of  r.  Then  by  the  perturbed 

'X/ 

graph  based  on  {v.| ,  v2»  w^  ,  w^)  ,  we  mean  the  graph  r  whose  vertices  are 
the  same  as  those  of  r,  and  whose  edge  set  is 

(Ej.  -  {(v1,v2),  (w-|  ,w2) } j  U  {(v1,w1),  (v2,w2)}  .  Here  Er  is  the  edge  set 

<\/ 

of  T.  We  say  also  that  T  is  the  result  of  a  perturbation  on  T.  These 

'b 

modifications  are  precisely  the  X-changes  defined  in  [T-S].  T  is  regular 

of  degree  d  if  r  is.  We  will  view  these  perturbations  as  "small" 

change  in  the  graph,  and  move  in  directions  which  improve  a  certain 

functional  which  we  will  define  below. 

Graphs  will  be  encoded  by  their  "incidence  matrices."  We  number  the 

vertices  of  the  graph  r,  {v^  ,  ...,  v^}.  By  the  incidence  matrix  I(r), 

we  mean  the  matrix  [a..],  where  a-.  =  1  if  v.v.  is  an  edge  of  T,  or 

1  J  •  J  1  J 

i  =  j,  and  a.  .  =  0  otherwise.  One  useful  property  of  I ( r )  is 
*  %) 

k 

The  (i,j)-th  entry  of  I(r)  is  the  number  of  paths  of  length  <  k 


from  v.  to  v.  in  T. 


Consequently,  we  have  that  the  diameter  of  T  is  the  least  value  of 
k  for  which  all  the  entries  of  I (r)  are  non-zero.  It  is  this  criterion 
which  is  used  in  the  algorithm  to  compute  the  diameter,  since  matrix  powers 
are  readily  computable  by  a  machine. 
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It  turns  out  that  the  diameter  alone  is  itself  not  a  sufficiently  sensi¬ 
tive  invariant  for  purposes  of  the  algorithm.  Specifically,  there  are  too 
many  graphs  for  which  no  perturbation  results  in  an  improvement  of  the 
diameter.  Consequently,  using  only  thp  diameter  as  a  functional  to  be 
optimized,  the  algorithm  is  frequently  unable  to  find  a  graph  with  even 
reasonable  diameter.  To  see  how  to  improve  matters,  we  need  a  definition. 

Let  TL  denote  the  integers.  We  wish  to  define  an  ordering  on  2n  ,  called 
the  lexicographic  ordering.  For  n  =  l,  the  lexicographic  ordering  on  Zn  =  TL 
is  just  the  usual  ordering  on  the  integers.  For  n  >  1 ,  we  suppose  the 
ordering  is  already  defined  for  all  m  <  n.  We  write  2Zn  -  TL  x  Zn-^  ,  and 
define  the  ordering  inductively  by 

z  <  z'  or 

(z,w)  <  ( z '  ,w 1 )  - 

z  =  z'  and  w  <  w'  . 

n  i 

for  z  e  TL  ,  w  €  TL  .  Now,  we  will  associate  to  every  graph  its 
"diameter  vector."  First,  for  a  positive  integer  H,  we  define  ot^r)  to 
be  the  number  of  zeros  in  the  £-th  power  of  I(r).  Of  course,  if  %  >  A(r), 
at(D  =  0.  The  diameter  vector  is  now  simply  the  vector  -fs~~rrew-simp-I^  tto — ■ 
-vector — 

^A(r)(-k) ,  ak-1Cr),  ak_2(r),  ....  a2(r)j  . 

We  will  denote  this  by  v(r).  We  order  these  vectors  as  follows: 

(Mr),  ak_l(r)*  a2(r))  “k'.i^)*  •••»  a2(r)) 

A(r)  <  A(r')  or 
A(T)  =  A(T')  and  a(r)  <  a(T') 


otgfr)  \  ,  and  the  ordering  is 
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Here,  a(r)  denotes  the  vector  ^a^^Cr),  .... 
the  lexicographic  ordering. 


The  algorithm  now  proceeds  as  follows.  We  say  a  4-tuple  (v^  ,  v2,  w1 ,  w2) 
is  T-admissible  if  v^v2  and  w^w2  are  edges  of  T,  and  and  v2w2 

are  not.  To  a  T-admissible  4-tuple  we  may  associate  a  perturbation  of  T, 


as  defined  above.  From  a  fixed  initial  graph  T,  T-admissible  4-tuple*?  arc 


generated,  and  the  associated  perturbations  are  applied.  This  continues  for 


a  large  number  of  steps,  until  the  initial  graph  is  presumed  randomized. 
The  4- tuple/  are  generated  using  a  random  number  generator.  The  fixed 
initial  graph  (in  the  trivalent  case)  is  an  n-cycle  with  antipodal  points 


connected.  After  this  is  done,  the  steps  are  as  follows: 


(I)  Select  a  T-admissible  4-tuple  at  random* 

'll  *\j 

(ii)  Compute  V(r),  where  r  is  the  graph  obtained  by  applying  the 
perturbation  associated  to  the  4- tuple  constructed  in  (I). 

If  v(T)  <  v(T) ,  set  T  =  T.  Repeat  step  (I). 


The  perturbations  are  selected  at  random,  since  it  was  found  that  a 
simple  ordering  of  perturbations  tended  to  bias  the  algorithm  toward 
particular  graphs. 

We  compare  our  algorithm  to  that  devised  in  [T-S].  Our  perturbations 
are  precisely  their  X-changes,  but  the  functional  we  optimize  is  much  more 
sensitive.  Theirs  consists  only  of  A(T)  and  a^_-|(T). 

We  now  summarize  the  results  of  the  application  of  our  algorithm. 

By  the  use  of  the  algorithm,  it  has  been  possible  to  improve  substantially 
most  of  the  densest  known  graphs.  We  give  our  improved  version  of  the  table 
constructed  in  [LFQSU].  d  denotes  the  degree  of  the  graph,  k  the  diameter. 
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The  (d,k)  entry  is  the  largest  known  graph  with  diameter  k  and  degree  d. 

Our  entry  is  listed  above,  the  parenthesized  value  below  is  the  value  from 

[LFSQU].  One  asterisk  indicates  that  the  graph  is  probably  optimal.  Two 

asterisks  indicates  that  it  is  obtained  from  [B]  using  the  result  cited  in  §1. 
The  results  obtained  from  this  algorithm  are  in  some  cases  surprising. 

We  list  some  of  the  qualitative  properties  observed. 

(a)  Many  graphs  obtained  by  random  generation  of  graphs  improved  values  in 
the  older  version  cf  the  table  in  [LFQSU],  given  by  Storwick  [S ].  This 
suggests  that  one  is  further  from  optimal  than  was  previously  thought. 

(b)  By  evaluating  the  eigenvalues  of  the  incidence  matrices  arrived  at  by 
the  algorithm,  it  was  found  that  there  are  many  distinct  "local  minina" 
(i.e.,  graphs  for  which  no  perturbation  improves  the  diameter  vector) 
for  the  diameter  vector.  This  contradicts  the  suggestion  made  in  [T-S], 
that  one  tends  to  arrive  at  a  global  optimum  from  all  starting  points. 

It  seems  that  the  algorithm  in  [T-S]  suffers  from  two  deficiencies. 

First,  their  objective  functional  for  minimization  is  not  sufficiently 
sensitive,  as  we  observed  above.  Secondly,  their  perturbations  are  done 
in  fixed  sequential  order,  which  severely  skews  their  results.  We  have 
overcome  this  difficulty  by  randomly  selecting  the  perturbations  at  each 
stage. 

(c)  Although  our  algorithm  is  efficient,  it  seems  that  substantially  larger 
networks  could  be  studied  if  our  diameter  routine  were  modified  to  use 
tne  so-called  "Dijkstra  algorithm,"  which  would  speed  up  the  diameter 
calculation  substantially. 


(d) 
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Although  the  algorithm  is  an  improvement  over  all  previous  heuristic 
algorithms  for  this  problem,  ’t  is  unable  to  find  many  known  dense 
graphs,  arising  from  systematic  constructions.  The  reason  for  this 
seems  to  be  the  "denseness"  of  the  set  of  local  optima  in  the  set  of 
ail  graphs  of  degree  d,  and  the  relative  sparseness  of  the  so-called 
vertex  transitive  graphs  therein.  It  seems,  therefore,  that  it  would 
be  desirable  to  design  and)  si  gorithm  which  operates  entirely  inside  a 
collection  of  vertex  transitive  graphs,  possibly  within  the  Cayley 
graphs  (see  ?D).  Using  the  Dijicstra  diameter  algorithm,  and  an 
efficient  description  cf  many  groups,  such  an  algorithm  should  be 
construct! bl - .  Moreover,,  it  'would  allow  much  larger  networks  to  be 
studied,  since  the  diameter  calculation  for  vertex  transitive  graphs 
is  substantially  shorter- than  that  for  arbitrary  graphs. 

A  copy  of  the  code  for  the  implementation  cf  this  algorithm  is  enclosed. 


as  well  as  documentation  fer  it. 
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( C )  Modification  or  the  Algorithm  for  More  Sensitive  Measures 

£ 

In  view  of  the  remarks  in  SI  which  identify  the  measures  A^,(r)  as  the 

p 

diameter  of  an  associated  grapn  one  can  in  principle  study  these 

£ 

measures  using  the  algorithm  discussed  in  SB.  However,  the  graphs  F^  are 
usually  too  large  for  this  procedure  to  be  practicable.  Cur  current 

implementation  of  the  algorithm  will  accept  only  graphs  with  fewer  than 

o  o 

10,000  points,  end  usually  is  larger  than  this.  For  the  measures  A.  , 

we  therefore  propose  the  use  of  a  “dual  algorithm"  based  on  girth,  which  is 
much  simpler  to  ccmoute  (see  §A). 

The  modified  girth  algorithm  is  identical  to  the  previous  algorithm, 

except  that  the  objective  Functional  is  altered.  For  a  graph  r,  we  define 

its  girth  vector  to  be  y(F)  =  (v.j(r),  v,(F),  ...,  vk(r),  ...  j.  This  is 

ordered  by  the  lexicographic  ordering,  ar.d  the  algorithm  proceeds  just  as 

before,  except  that  we  r.cw  accept  a  perturbation  if  ~.t  increases  y(f). 

0 

Applying  this  algorithm  to  F^  should  produce  heuristic  results  which 
improve  these  measures. 


1b 


(D)  Vertex-Transf ti ve  Graphs 

One  desirable  feature  of  a  graph  to  be  used  as  a  processor  array  is 
that  the  description  be  as  simple  and  compact  as  possible.  Thus,  the  graphs 
produced  by  a  heuristic  al gorithm  will  in  general  not  be  satisfactory  from 
this  point  of  view,  and  for  this  purpose  it  would  be  useful  to  restrict 
oneself  to  a  class  of  graphs  having  a  compact  description. 

One  such  family  is  the  collection  of  so-called  Cayley  graphs.  Let  us 
define  these.  Let  G  be  a  finite  group,  ana  let  £2  £  2  be  a  subset, 
closed  under  inversion,  so  that  £2  =  £2~^ .  Then  v/e  define  the  Cayley  graph 
associated  to  the  pair  (G,£2),  r(G,£2),  to  be  the  graph  whose  vertices  are 
the  elements  of  G,  and  whose  edges  are  pairs  (e,ug)  for  <  £2.  As  an 
example,  if  G  =  ZZ/n,  the  cyclic  group  with  n  elements,  and  P.  =  {T,T-^}. 
where  T  is  a  generator  of  G,  the  associated  graph  is  tfie  cyclic  graph  of 
size  n.  A  useful  property  of  r(G,P)  is  that  its  automorphism  group  is 
transitive  on  the  vertices,  i.e„ ,  for  every  pair  {x,y5.  of  vertices  of 
r(G,£2),  there  is  an  automorphism  a  of  r^G.P)  sc  that  a(x)  =  y.  This 
is  clearly  the  case  sines  the  right  G-acticn  on  C  provides  an  action  of 
G  on  the  graph,  which  is  clearly  transitive  on  the  vortex  set.  This  is 
a  useful  property,  since  it  means  that  the  diameter  may  be  computed  by 
finding  the  ooints  of  maximal  distance  from  one  given  point.  Also,  routing 
algorithms  for  these  networks  are  compactly  described,  since  one  must  only 
find  optimal  paths  sorting  at  one  given  point. 

One  important  proposed  architecture,  the  cube-connected  cycles  of  [P- V 

are  in  fact  of  this  roiv.  For,  consider  the  semi-direct  product 

G  =  2Z/n  x  2Z/2n  ,  p{T){x-»  x  }  =  (xp,  ^ ,  ....  x^).  Thus  G  has 

P 


r  *i 
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clementc  (m,v),  where  m  <=  Z/n  ,  and  v  C  (Z/2)n  ,  and 
(m,v)(m’,v')  =  (m+m1 ,  p(m')v  +  v‘).  It  is  an  easy  calculation  to  see  that 
if  ft  =  {(1,0),  (-1,0),  (Oje^)},  where  e^  =  (1 ,0,0, . . .  ,0) ,  then  r(G,ft)  is 
in  fact  isomorphic  to  the  graph  associated  with  the  cube-connected  cycles. 

The  diameter  of  the  cube-connected  cycles  is  known  to  grow  as  5/2  log,,  (K), 
where  K  is  the  number  of  vertices  in  the  graph. 

We  should  remark  here  that  large  girth  (and  hence  small  diameter)  in 
Cayley  graphs  is  associated  with  non-commutativity  of  the  group  in  question. 
This  being  the  case,  the  simple  groups  seem  to  be  natural  candidates  to 
produce  efficient  graphs.  To  see  that  this  is  in  fact  the  case,  we  study 
the  diameters  of  r(G,ft)  for  certain  choices  of  ft,  and  G  =  S^,  and  observe 
that  for  large  n,  they  approximate  the  Moore  bound.  The  order  of  the 
symmetric  group  Srt  is  n'.  .  Let  ft  c  $  be  the  set  of  all  cycles  of 
length  <  k.  We  propose  to  compare  the  diameter  cf  r($n ,ft,. )  with  its 
associated  Moore  bound.  First  we  observe  that  ]ftv|  -  (l)}  +  2(^)  +  6(^)  +  ... 

+  (k-1 )  1  (j^) ,  by  a  simple  counting  argument.  Thus,  the  degree  of  r(Sn>ftj,) 


k 

is  d  *-  Z  (j-1)  1  (”) .  Thus  for  large  n,  the  Moore  bound  for  7[7  .ft,.' 
0=2  J  n  K 


is  logcl_1(nl)  =  ^fferr  .  By  Stirling’s  formula,  iim  =  1,  so  fe 

'  1  n 


n  ?,n  n 


large  n,  the  Moore  bound  is  approximated  by  \-,-y  .  But,  aaa.i,  fov 

.oH  \  U-  \  / 


/  k 

large  n,  £n(d-l)  =  Jin  (  I  (j-1) 

V"2 


is  approximated  by  S.n %  k  ?.n  n. 


Her.ce,  for  large  n,  the  Moore  bound  for  r(Sn,ftj,)  is  approximated  by  u/.:. 
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The  diameter  of  (S^,^),  on  the  other  hand,  tends  to  r./k-l,  as  one 

readily  computes  in  S  :  Consequently,  the  diameter  of  (S  ,ft.)  is  within 

n  n  k 

a  factor  of  k/k-l=  1  +  of  the  Moore  bound.  As  k  becomes  larger, 

we  are  able  to  approximate  equality  arbitrarily  closely.  This  shows  that  S 

is  a  plausible  candidate  for  further  study. 

We  will  now  show  how  to  modify  the  cube-connected  cycles,  using  group 

theoretic  methods,  tc  provide  on -infinite  family  of  vertex  transitive 

trivalent  graphs,  whose  diameter  is  substantially  smaller,  but  which  have 

all  the  desirable  regularity  properties  of  the  cube-connected  cycles. 

Let  5  =  Z/n  x  (Z/2)r',  as  before.  Note  that  G  contains  a  central 
n  p  .  r. 

element,  namely  the  vector  (0,(1 ,1 ,... ,1 )}.  Following  our  intuition  con¬ 
cerning  the  relationship  between  non-commutativity  me  small  diameter,  we 

eliminate  the  central  element  by  simply  factoring  it  out,  Call  the 
n,  o, 

quotient,  group  and  let  Q  denote  the  image-  of  H.  Then  we  claim 
that  the  diameter  cf  ?(Gn,w)  grows  as  2  logoi'),  where  K  is  the  number 
of  vertices,  an  improvement  over  the  -diameter  5/2  log«,(K)  obtained  for 
the  cube-connected  cycles. 


'"0  V  ^ 

Proposition  The  diameter  of  r ( G_ -JT2)  prows  as  2  loq0(K) 

_ _ M _  (,. 

_/w  0/ 

Proof  By  taking  inverses,  we  can  clearly  consider  the  graph  r f Gn ,Q ) . 
where  ( a , g co)  is  an  edge  for  n  t  fi.  An  element  of  is  giver,  by  an 
ordered  pair  (m,v),  m  k  Z/.i  ,  v  6  V  (Z./C)n/(1 , . . .  ,1 } .  The  multiplica- 
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the  first  coordinate  of  a  vector  v  S  V  -  (Z/2)n.  For  any  v  C  V,  one  may 

f\j 

lift  v  to  an  element  v  of  V,  so  that  the  number  of  non-zerc  coordinates 
in  v  is  <  n/2,  for  if  one  lift  v  does  not  have  this  property,  then 
v  +  (1,1,...,  1)  does.  Given  (m,v),  select  v  as  above. 

The  algorithm  now  proceeds  as  follows. 

(I)  Initialize  a  counter  a  at  n-1 . 

(II)  Is  x,  ( v)  =  0  ?  If  yes,  proceed  to  (IV),  if  no,  proceed  to  (111). 

(Ill)  Multiply  by  e^.  Proceed  to  (IV). 

(IV)  Multiply  by  T.  oi  -  1  -►  a  .  Proceed  to  (V). 

(V)  Is  c:  =  0  ?  If  so,  proceed  to  (VI).  If  net,  return  to  (IT). 

(VI)  If  m  r  0,  multiply  by  Tq,  where  q  -is  the  number  of  minimal 

absolute  value  congruent  to  -m  mod  n.  i-.'otc  that  < q (  <  n/2.  Quit. 

Since  v  has  at  most  n/2  1‘s,  we  only  multiply  by  e.  a..  most  n/2 
times.  Thus,  the  total  number  of  steps  is  at  most  (n-1)  +  n/2  -  n/2  -  2n  -1. 
For  odd  n,  it  is  at  roost  2n  -  3,  which  is  of  the  same  order  as 
2  log,(n  2n_1). 

^  1.1 


If  one  forms  the  quotient  of  this  graph  by  the  equivalence  relation 
(r.:,v)  (m1  ,v) ,  V  m,rri'  s  71/ n  ,  one  obtains  a  non-regular  family  of  graphs 

whose  diameter  grows  as  3/2  1092(0),  which  is  comparable  to  that  obtained 
in  [L-S'J,  and  for  which  the  routing  algorithm  is  much  simpler.  Finally, 


an  alternative  version  of  this  construction  is  given  by  forming  the 
(n-1  )-cube,  inserting  n-cycles  at  every  vertex  so  that:  the-  incoming  edges 


etch  connect  at  distinct  vert'ces,  and  connect 


i.i'.e  *eiivnni ng  vertex  to  the 
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(E)  A  Layout  for  the  Modified  Cube-Connected  Cycles 

•  In  the  paper  [J>-V],  two  layouts  are  proposed  for  the  cube  connected 

cycles,  one  slightly  more  efficient  than  the  other.  By  combining  these  two 

layouts,  we  obtain  a  layout  for  the  modified  cube  connected  cycles.  The 

area  of  the  layout  grows  as  2/2  times  the  area  of  the  cube  connected  cycles 

with  the  same  number  of  nodes,  and  has  communication  time  roughly  4/5  times 

that  of  the  cube  connected  cycles.  We  give  the  layout  for  the  case  n  =  5. 

4 

corresponding  to  5*2  =80  nodes.  It  will  be  clear  from  the  diagram  how 
to  extend  to  the  general  case. 
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makefile 


makefile 


CFLAG5 =  -0 
search  : 

search.o  : 

zeroes.o  : 

eigen.o  : 
geneval.o: 


search.o  zeroes.o  eigen.o  geneval.o 

cc  S(CFLAGS)  search.o  zeroes.o  eigen.o  geneval.o  — lm 

declar.h  search.c 
cc  S(CFLaGS)  — c  search.c 

declar.h  zeroes. c 
cc  S(CFLAGS)  -c  zeroes. c 

declar.h  eigen.c 
cc  S(CFLAGS)  — c  eigen.c 

declar.h  geneval.c 
cc  S(CFLAGS)  -c  geneval.c 
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declar.h 


^(define  MOD  % 

^define  Intsize  (8  *  sizeof(long)) 

^define  Maxsize  256 
#define  Maxreduced  Maxsize /(Intsize) 
jjldefine  YES  1 
^define  NO  0 

^define  SCREENSZ  80  /*  the  width  of  the  screen;  for  I/O  control  •/ 
^define  EV3IZE  6  /•  the  width  of  the  field  of  an  e.v  display  •/ 

#define  DISPCT  (SCREENSZ  /  EVSIZE) 

extern  int  (•multiply)(); 
extern  int  raultgraph(); 
extern  int  multlgraph(); 
extern  int  mult2graph(); 

extern  int  (•ne?-tgraph)(); 
extern  int  se  ntial(); 

extern  int  rr  _move(); 

extern  int  (*svaluategraph)(); 
extern  int  eval_diamvect(); 
extern  int  eval_passa!l(); 

typedef  long  Graph[Maxsize][Maxreduced]; 
typedef  double  Matrix[Maxsize][Maxsize]; 
typedef  double  Vector[MaxsizeJ; 
typedef  long  Passgraph[][Maxreduced]; 

extern  int  numpts;  /•  the  number  cf  pts  in  the  graph  •/ 

extern  int  degree; 

extern'  long  place[Maxsize]; 

eXtem  long  graphA[Maxsize][Maxreduced1: 

extern' long  graph[2][Maxsize][Maxreducedj; 

eXtem  double  matrix[MaxsizeJ[Maxsize]; 

extern  double  ewec[Maxsize]; 

eXtem  int  vectorA[Maxsize],  vectorB[Maxsize]; 

extern  int  *vectA,  *vectB,  «temp; 

extern;  int  diamA.  diamB; 

extern  int  xl,  x2,  x3,  x4;  /*  coords,  for  the  current  perturbation  •/ 

extern  struct  Status  f 

int  starting;  /*  =YES  =>  start  up  procedures  for  the  program  */ 

Int  coctnue;  /•  =N0  =>  don't  look  for  any  more  graphs  */ 

int  newgraph;  /•  =YES  =>  last  perturbation  was  accepted  •/ 

long  graphcount;  /•  4  of  new  graphs  produced  •/ 

hit  printing;  /•  =NC  =>  no  printing  of  intermediate  results  */ 

Int  diamwatch;  /•  to  watch  diam.  results  for  unsuccessful 

comparisons  set  diamwatch  =  YES  •/ 
int  backtracking;  /•  =YES  =>  backtracking  in  effect  •/ 
int  eigenvalues;  /•  =YES  =>  find  all  eigenvalues  */ 

1  status; 

eXtemdcuble  Sqrt2mv; 

extern  long  ncxtprint,  increment,  incincrement,  incfactcr; 
extern  long  lowerlimit,  uppcrlimit; 
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/include  <math.h> 

/include  "declar  h" 

double  matrix[Maxsize][Maxs}ze]; 
double  evvec[Maxsize]; 

compar(numl,  num2)  COTTip'LT 

double  *numl,  *num2; 

•/♦  this  procedure  compares  two  real.. numbers,  returning  -  I  if  numl  > 
nurr.2,  0  if  num  1  =num2,  and  + 1  if  numl  <  nurr.2.  Big  deal.  (It  is  needed 
.  for  qsort  in  cyclic  jjacobi.)  */ 

if  (*numl  >  “nurci2) 
retum(-l);  " 

else  If  (•numl  =■=  *num2) 
i-eturn(C); 
else 

return(l); 

i 


. . . . . 

cyclic-^acobi(graphl)  Cyclic CLCOb'L 

Passgraph  graphl; 

•/♦  this  procedure  tak.es  graph 1,  which  is  the  "augmented" incidence  matrix 

■  of  a  graph,  strips  crxi  the  diagonal  l's,  leaving  the  usual  incidence  matrix, 

:  transforms  the  "bit  level" matrix  to  one  with  real  entries  (only  the  ripper 

.  half  of  the  matrix  is  used  anywhere),  and  then  uses 

•  the  cyclic  Jacobi  method  to  find  the  eigenvalues,  ft  also  comput  20171*3 

•  measure  of  the  degeneracy  (still  to  be  determined  at  this  printing). 

"cvvecl"ic  the  quicksorted  diagonal  of  the  matrix,  ft 

■  seems  very  likely  that  for  large  graphs  better  methods  using  the  sjur-sify 
=  of  the  incidence  matrix  would  be  valuable.  */ 

l 

/define  EPSILON  0.000001 
£  define  HER3D  W0 
double  error; 
int  itnum; 
double  *ai; 

Ions  •£'! 
int  i,  j; 

lor(i=0;  i  <  numpls,  i++) 

t 

ai  =  matrix[i]; 

•(ai  +  i)  =  C.C; 
gi  =  graph  l[ij; 
fer(j=i+l;  j  <  numpts;  j-*-+) 
if  (*(g>  +  (j/Intsize))  &  p(ace(j]) 

*(ai  +  j)  =  1.0; 
else 

•(ai  i-  j)  -  0.0; 

I 

error  =  1.0; 

for  (jtnum  ~  0;  error  >  EPSILON;  itnum-*-+) 

t 

cytaraticr.f  matrix,  deerror); 
if  (itnum  >  ITERBD) 

4 

t 

prir.tfC'CJ  Alg.  looping  too  many  times 
return, 

} 
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l  oyi:  i  of  cIqh.v  c 


eigen.c 


i 

for(i=0;  i  <  numpts;  i-“~- ) 
ewec[i]  =  matnx[i][i]; 

qsort(ewec,  numpts,  irixeol(dauble),  compar); 
evprint(ewec); 

t 

. . * . . . 4.... ....4.. . . 

cji*eration(matnxl,errp.t) 

Matrix  matrix  1; 
double  *errpt; 

/•this  procedure  performs  one  iteration  of  the.  cyclic  -Jacobs  algorithm 
on  the  matrix  matnx2.  The  error  is  the  sum  of  the  squares  of  the 
off-diagonal  elements.  .Voticc  that  only  the  upper  half  of  the  matrix  is 
ever  used.*/ 

I 

f  define  DELTA  O.OOCOOaOOOOO  1 
int  i,  j,  k; 
double  q,  v; 
double  alpha,  beta; 
double  alphasqr; 
double  aij,  aijsqr; 
double  all.  ajj; 
double  dmi,  dmj; 
double  temp; 
double  *ai.  *aj; 

for(i=0;  i  <  numpts; 

i 

ai  =  matrix1[i]; 
for(j=i-*-l;  j  <  liumeta;  j-r-r) 

l 

aij  =  *(ai  j); 
aijsqr  =  pcvr(aij.  2.0); 
il  (aijsqr  >  DELTA) 

l 

aj  -  matrix  l[;j; 
aii  =  *(ai  +  i); 
ajj  =  .(aj  +  j); 
q  =  aii  —  ajj; 

v  =  sqrt(  (4.0  •  aijsqr)  +  pow(q,  2.0)  ); 

?r  (q  «=  c.o) 

beta  =  alpha  -'SqrtSinv; 
alphasqr  =  0.5; 

efcie 

{ 

alphasqr  =  tab$(v  +  q)  /’(2.0  *  v); 
alpha  =  sqrt(alphasqr); 
beta  =  aii  /  (v  *  alpha); 

V 

fer(k=0;  k  <  i;  k++) 

dmi  =  matrix  l[k][ij; 
nmj  -  matrixt[kj[jj; 

metrixl[kj[i]  —  alpha  •  dmi  beta  •  dmj; 
vnatri;:!  [!;)[;]  -  alpha  *  d;nj-  beta  *  dmi; 

f  , 

1  ';mp  =  (q'  alphasqr)  +  (3.0  *  nirhe  *  beta  •  ai};; 

“'(ai  -*•  i)  =  ajj  J-  temp; 


eigen. c 

...cyclic— jacobi 


cjiteration 
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eigen,  c 


eigen. c 


for(k=i+l;  k  <  j;  k — ) 

dmi  -  »(ai  +  k); 
dmj  =  matrix  l  [k][jj; 

•(ai  +  k)  =  alpha  •  ami  +  beta  •  dmj; 
roatrixl[kj|.j]  =  alpha  •  dmj  -  beta  •  dmi; 

t 

♦(ai  +  j)  =  0.0; 

•(aj  +  j)  =  aii  -  temp; 
fcr(k=j+l;  k  <  numots;  k++) 

{ 

dmi  =  *(ai  +  k); 
dmj  =  »(aj  +  k); 

•fai  +  k)  =  alpha  ♦  dmi  -  beta  •  dmj; 

*(aj  -r  k)  =  alpha  •  dmj  —  beta  *  dmi; 

\ 


•errpt  =  0.0; 

for(i=0;  i  <  numpts;  i -»■-*-) 
for(j=i+i;  j  <  numpts;  j+-r) 

♦errpt  +=  pow(niatrij;l[i][j],  2.0); 

1  .  . 

. . . . 


evprint(ewec2) 

Vector 

/♦  this  procedure  prints  r he  eigenvalue  information  obtained  by 
cyclic _jacobi,  above.  */ 

A 

int  i; 

printf("Number  of  r.odes:%ld,  Degree:  TSld,  ’’.numpts.  degree); 
nrintfC’Eigenvaiues  of  the  incidence  rnatrix:\n7); 
for(i=0;  i< numpts;  i++j 

prin.tf('’%7.3f%s,,,ewec2[i]1((i  MOD  10  ==  9)jj(i==numpts  -  1))  ?  "\c"  : 


...cj  iteration 


evvrin 

j. 
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SO 


geneval.c 


geueval.c 


^include  "declar-h" 

generategraphQ  yG7lSTQ.LGCf7  CLph. 


if  (statusistarting  ==  YES) 

startgraph(graphA,  &numpts); 

return; 

(*nextgraph)(graphA); 


sequenliaKgraphi) 

Passgraph  graph  1; 

/•  takes  graph  l  and  identifies  the  next  sequential  perturbation, 

V  erf  or  ms  it,  and  sets  x) .  x2,  x3.andx4  to  identify  it  ♦/ 

i 

static  inti; 
state:  -mi  yS,  y4  ; 
static  rat  zl; 
static  iat  firsttinie; 

Adeline  13  ar.k  limit  30 

static  int  a  l[Backlixnit],  a2[Backlirnit],  aSfBackiimitj,  a4[BacklimitJ; 
static  int  c3[Back!in?it],  c4[Baeklim:ti.  j[Backlimit],  bl[Backlimit]; 
static  int  level  =  0; 


if  (status:r.er?£raph  --  YES) 

l 

firsttime  =  ^ES; 
if  (statas.graphcounl  !=  1L) 

\ 


al 

a2 

a? 

a4 

c3 

c4 


levei 

level 

level 

level 


=  xi; 
=  *2; 
=  x3; 
=:  x4; 


level)  =  y3; 


;cve 


*j  j 

1]  =  y4; 


j[leve(]  -  i; 
bl  [level]  =  zl; 

level  =  (leyel  -r  i)  JJCD^ar.klimit; 
goto  rezumi; 

» 

{ 

else 

goto  resume; 
level  -  0; 
for(  ; ; ) 

1 

for(xl=l;  xl  <  numots  -  1;  ?;!-+) 

* 

if  (firaUixnc  YES) 

$ 

firsttime  ~  nG; 
rl  =  xl; 

f 

r.loc  if  (xl 1 ; 

if  (status  !'3s*.:trr.ckJng  -----  V iif>) 

9 

f 

jn'intfCV/jr  (pc:isibly;  kwe/itu  gi 


sequv.ni.ial 
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printgraph(graphl,:numpts); 
cychc_jacobi(graphl); 
if  (diamA  ==  4) 

takepo(graphl,  vectA,  numpts,  <kdiamA,  numpts); 
status. graphcount — ; 
nextprint — ; 
mere  me  at  =  1L; 
incfactor-=  1L; 
incincrement  =  1L; 
lowerlimit  =  OL; 
upperlimit  =  10000G0L; 
diaminfo(vectA,  diamA); 

level  =  (level  —  1  ■+•  Backlimit)(MOD!Backlimit; 
xl  =  al  level]; 

x2  =  a2  level  ;  — 

x3  =  a3[level  ; 

X4  =  a4[level  ; 
y3  =  c3  level ; 
y4  =  c4[levelj; 
i  =  j[level]; 
zl  =  blflevel]; 

perturb(graphl,  xl.  x4.  x3,  x£); 

takepowers(graphl.  vectA.  numpts,  iediamA,  numpts); 
for  (i=0;  i  <  Maxsize;  i++) 

•(vectB  +  i)  =  •(vectA  +  i); 
diamB  =  diamA; 
goto  resume; 

I 

else 

J 

status.contnue  =  NO; 

return; 

i 

for  (x2=0;  x2  <  xl;  x2++) 
if  (graphl[xl][x2/(Intsize)]  St  placefx2]) 
for(y3=x'l  +  1;  y3  <  numpts;  y3+4) 
for  (y4=0;  y4  <  y3;  y4++) 
if  ((graphl[y3][y4/(Intsize)]  £qplace[y4]) 

&&;(xl  !=  y4):ScSt  (x2  !=iy4)) 
for(i=0;ii  <=  1;  i++) 

1  if  Ki==0) 

f 

x3  =  j3; 
x4  =y4; 

else 

x3  =;y4; 
x4  =y3; 

I 

if  (  !(graphl[x2][x3/(Intsize)]  it  place[x3]) 

'■Stic  !(graphl[xl][x4/(Intsize)]  St  place[x4])) 

perturb(graphl,  xl,  «2,  x3,  x4); 

return; 

l 


...sequential 


resume 


resuml 
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geneval.c 

...sequential 


l 


./■ 


/ 


mdm— move(graphl)  mdvn—m.OVe 

Passgraph  graph  1, 

/*  randomly  picks  the  next  perturb  alien,  performs  it,  and  records 
it  usina  xl,  x2,  x3,  x4  */ 

* 

inti,  j; 
for  (  ;  ;  ) 

i 

xl  =  randQ  MOD  nuinpts; 
j  ”  randQ  MOD  degree; 
for  (i=0*  i  <  numpts;  i++) 
if  (  (j  !=  xl)  kk  (graphl[xi][i/Intsizej  k  p!ace[i])  ) 
it  0**0) 

x2  =  i; 

break; 

? 

j 

else 

j — ; 

for  (x3=rand()  MOD  numpts;  (x3  =  xl)  ||  (x3  ==  x2);  x3=rand()  MOD  numpts); 
j  =  randQ  MOD  degree; 
fer  (i=C;  i  <  numpts;  i++) 
i I  (  (i  !=  >:3)  kk  (graphl[x3][i/Intsize]  <5c  place[i])  ) 

if  0  ==  0) 

i 

x4  =  i; 

break: 

p  . 

i 

else 

j — ; 

if  (  (x4  ’=  xl)  kk  (x4  !  =  *2)  kk  (’(graphl[x2][x3/Intsize]  k  place[x3j)) 
kk  (!(graDhl{xl][x4/Intsize]  k  Fiace[x4]))  ) 

1 

perturb(graphl,  xl,  x2,  x3,ac4); 

return: 

{ 

.else 

continue; 


. . *****.«**•**«<,•**.*•**•**•*•.•*••**•*«««****«*•«••*••••****••«*•#/ 

evel-diemvectQ  QVCll—d'LG.TTlV&ct 

\ 

inf  i; 

int  improved; 

static  Lnt  diamcompare; 

if  (status.starting  =-  1X5) 

i 

if  (numpts  <-  tntsiz**) 
multiply  =  multigraph; 
else  if  (numpts  <=  ,?.  ♦  Intsize) 
multiply  =  mullSgraph; 

else 

multiply  —  mullgraoh; 
zcroeounKgrapha,  riurr.pta, 
vectorA[l]  ~  i; 
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diamcompara  —  numpts; 

lakepowers(grapl:A,  vectA.  rumpti  &dian>A,  Jinmcnmpare); 
•(voctB  +  i)  =  i(vectA  1); 
diamB  -  dic'.rnA; 

i 

else 

i 

if  (status.diamwatch  ==  NO) 
diamcompare  =  diamA; 

takepowers(graphA,  vectB,  numpts,  ckdiamB,  diamcompare); 
compare(vectA,  vectB,  diaruA,  diamB,  ^improved); 
if  (improved) 

\ 

status. nevrgruph  =  YES; 
temp  ~  vectA; 
vectA  -  vectB; 
vectB  -  temp; 
diamA  =  diamB; 
if  (status. eigenvalues  =—  YES) 
cyclic_jacobi(graphA); 

else 

I 

pertu:rb(graphA,  :<1,  x4,  x3,  x2); 
status.newgraph  =  NO; 

> 

* 

\ 

if  (status.newgraph  ==  YES) 

I 

status.graphcount-r+; 

printccnt.roi(); 


/' 


eval— passallO 

/*  this  'considers  all  graphs  to  he  improvements  */ 

i 

status,  newgraph  =  YES; 

takepowers(graphA,  vectB,  numpts,  IrdiamB,  numpts); 
temp  =  vectA; 
vectA  =  vectB; 
vectB  =  temp; 
diairtA  =  diamB; 
if  (status.eigenvalues  ==  YES) 
cyc!ic_jacobi(graphA); 
status.newgraph  =  YES; 
status.graphcount-r-r; 
printcontrolQ; 


/•*•**«• 


perturb(graphl,  xxl,  xx2,  xx3,  xx4) 

Passgraph  graph  1; 
fntxxl,  xx2,  xx3,  xx4; 

/•  this  performs  the  perturbation 

(xzl;xx2),  (xx3,  xx 4)  ->  (xx2,  xxCI),  (xxl,  xx4) 


graphlfxxt 
graph  ii  xx2 
graph  llxx3 
graph  l[xx4 


xx2/(Ints’ze) 

xxl/(Intsize) 

xx4/(lntsizeV 

xx3/(Int3ize) 


&= 

U~ 

<fc= 

&= 


"■“(place 

'■"(place 

~(place 

~(plaee 


’xx2 

xxl 

xx4 

’xx3 


I); 

!): 


♦/ 


/ 


geneva^.c 

...cvaL—dic.mvect 


sval-passall 


; perturb 
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geneval.c 


geneval.c 

...perturb 


graphl[xxl][xx4 /(Intsize)]  |=  place[xx4' 
graphl[xx4]Lxxl /(Intsize)]  |=  place[xxl 
graphl[xx2]rxx3/(Intsize)]  |=  place  xx3 
graphl[xx3][xx2 /(Intsize)]  |=  place[xx2^ 


. . . . ••••••• . . 

compare(dumvecl,  dumvec2,  dumdiaml,  dumdiam2,  dumimp)  COTTipCLTR 

/•  This  procedure  compares  ( starting  at  the  bottom)  dummy  vecl  with 
dummy  vec2  and  if  2  is  better,  dumimp  is  set  to  1.  */ 
int  •dumimp; 

int  dumvecl[],  dumvec2[]; 
int  dumdiaml,  dumdiam2; 

l 

int  i  ; 

if  (dumdiam2  <  dumdiaml) 

f 

•dumimp  =  1; 

return; 

I 

•dumimp  =  0; 

for (i=dumdiaml;  i  >=  2;  i — ) 

l 

if  (dumvec2[i]  <  dumvecl[i]) 

•dumimp  =1; 

return;  . 

I 

else  if  (dumveel[i]  <  dumvec2[i]) 

•dumimp  =  0; 

return; 
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^include  <math.h.> 


^include  "declar.h" 

int  (*multiply)(!); 

Int  multgraph(); 
int  multlgraph(); 
int  mult2graph{); 

int  (*nextgrapH)(); 
int  sequential); 
int  rndm_moveX); 

int  (*evaluategraph)(); 
int  eval_diamvect(); 

int  eral_paasaiK); 


intnumpts;  S+  the  number  of  pis  in  the  graph  •/ 

int  degree; 

long  place[Maxsize]; 

long  graph[2][Maxsize][Maxreduced]; 

long  graphA[Maxsize  ][Maxreduced]; 

int  vectorA[Maxsize],  vectorB[Maxsize]; 

ml  *veclA,  •vecTB.  *temp; 

int  diamA,  diamB; 

int  xl.  x2.  x3,  x!4;  /•  coords,  for  the  current  perturbation  •/ 

struct  Status  status  =  J 
YES,  YES.  YES.  OL/YES.  NO.  YES.  NOJ; 

/• 


int  starting; 
int  contnue; 
int  newgraph; 
long  grapheount; 
int  printing; 
int  diamwatch; 

int  backtracking; 
.  int  eigenvalues; 

•  / 


/  ~  YES  -  >  start  up  procedures  for  the  program  / 
/  -=NO  =>  -don't  -Look  for -any  more -graphs  / 

/  -YES  =>  last  perturbation  was  accepted  / 

(/  /  of  new  graphs  produced  / 

/  =NO  =>  no  printing  of  intermediate  results  / 

V  to  watch  diam.  results  for  unsuccessful 
comparisons  set  diamwatch  =  YES  / 

V  =  YES =>  backtracking  in  effect  / 

/  -  YES  =>  find  all  eigenvalues  / 


doable  Sqrt2inv; 

long  nextprint  =  1L.  increment  =  1L.  incincrement  =  OL,  inefactor  =  lL; 
long  lowerlixnit  =  1L,  upperlimit  =  1000000L; 


mainO 

i 

inti; 

evaluategraph  =  eval_diamvect; 
nextgraph  =  sequential; 

srand(l); 

Sqrt2inv  =  1.0  /  sqrt(2.'0); 
vectA  =  vectorA; 
vectB  =  vector  B; 

place[0]  =  1L; 
for(i=l;  i  <  Iiitsize;  i++) 
placeft]  =  2“  place^i— 1]; 
for  (i=Intsize;  i  <  Maxsize;  i++) 
place[i]  =  place[i  HOD  (Intsize)]; 

for  (  ,  status. contnue  ==  YES; ) 
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i 

generategraphQ; 
if  (status. contnue  ==  NO) 

break; 

(*evaluategraph)(); 

printf("Our  favorite  graph:\n’’); 
printgraph(graphA,  numpts); 
printf("\n''); 

takepowers(graphA,  vectA,  numpts,  &diamA,  numpts); 
diaminfo(vectA,  diamA); 

I 

. . . . *..../ 

printcontrol() 

I 

int  reply,  okreply; 

if  (status. printing  ==  NO) 

return; 

if  (status.graphcount  ==  nextprint) 

if  (status.graphcount  >=  lowerlimit) 
diaminfo(vectA,  diamA); 
nextprint  +=  increment; 

increment  =  incfactor  •  increment  +  incincrement; 

I 

if  (status. starting  ==  YES  ||  status.graphcount  >=  ujjperlimit) 

i 

status.starting  =  NO; 

for  (okreply  =  NO;  okreply  ==  NO; ) 

printf("How  would  you  like  the  intermediate  resullts  printeid?\n''): 

printf("No  intermediate  results  printed  (0).\n"); 

printf("Doubling  of  the  interval  between  prints  (().\nf); 

printf(”A  print  at  every  change  of  diameter  (2).\iT); 

printf("A  semi-custom  print  interval  (3).\n"); 

prlntf("Print  at  each  rtep(4).\ri"); 

printf("VThat  is  your  choice?  ;  "); 

scanf(''%d",  Areply); 

printf("\n'’); 

switch  (reply) 

i 

case  0: 

status  .printing  =  NO; 
nextprint  =  0L; 

return; 
case  1: 

nextgraph  =  sequential; 
evaluategraph  =  eval_diamvect; 
nextprint  =  status,  grephcourit  +  1; 
lowerlimit  =  0L; 
upperlimit  =  1000000L; 

Increment  =  1L; 
incincrement  =  0L; 
incfactor  =  2L; 
okreply  =  YES; 
break, 
case  2 

printf("option  2  is  not  available  yet\n”); 

break, 
case  3. 

printf("The  current  graph  count  is:  5£ld\n'\  status. graphfcbunt); 


printcontrol 
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...print  control 

printf(”Enter  the  count  for  the  start  of  printing:  "); 
scanf("751d",  &lowerlimit); 
nextprint  =  lowerlimit; 

printf("Enter  the  count  for  reconsidering  print  control  info:  "); 
scanf("%ld",  Aupperlimit); 

printf("After  printing  starts,  printing  will  occur  at\n”); 

printfC'increments  of  'a'. to  start  with,  and  changed  accardingNn”); 

printf("to  the  falling  rule:  new  =  'b'  •  old  ■+■  'c' .  ”); 

printf("at  each  step.Nn”); 

printf("Please  enter  a,  b,  c  an  this  line:  "); 

scanf("%ld  %ld  7.1d",  <kincrement,  icincfactor,  &incincrement); 

printf("Which  control  mechanism  would  you  like  for  "); 

printf("generatmg  new  graphs?\n"); 

printf("(l)  sequential\n''); 

printf("(2)  sequential  with  backtracking^"); 

printf("(3)  random\n"); 

printf("Yaur  chaice:'"); 

scanf("%d",  &reply); 

switch,  (reply) 

f 

case  1: 

nextgraph  =  sequential; 
status.backtracking  =  NO; 

break; 
case  2: 

nextgraph  =  sequential; 
status. backtracking  =  YES; 

break; 
case  3: 

nextgraph  =  nndm_mave; 

break; 

default: 

printf("Bad  aptiamselectedNn"); 

break; 

printf("Which  evaluation  function1  should  be  used  far  accepting"); 

printf("new  graphs?\n"); 

printf("(l)  use  Che  zeracounl  vector\n"); 

printf("(2)  accept  all  new  graphs  Nn"); 

printf("Yaur  choice;”'); 

scanf("75d",  <5creply); 

switch  (reply) 

l 

case  1: 

evaluategraph  =  eval_diamvect; 

break; 
case  2: 

evaluategraph  =  eval_pa3sall; 

break; 

default: 

prlntf("Bad  aptiamselectedNn"); 

break; 

l 

printf("Indicate(  whether  you 'are  happy  (1)  or  unhappy(2)  with  "); 

printf("your  choices:  "); 

scanf("%d",  Areply); 

okreply  =  (reply  ==  1)  ?  YES  :  NO; 

break, 
case  4; 

nextgraph  =  sequential; 
evaluategraph  =  eval_diamvect; 
nextprint  =  stalus.graphcount  +  1; 
lowerlimit  =  OL; 
upperlimlt  =  1000000L; 
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increment  =  1L; 
incincrement  =  OL; 
incfactor  =  1L; 
okreply  =  YES; 

break; 

default: 

printf("why  did  you  type  in  %d  —  please  start  aver";  reply); 

break; 


printgraph(graphl,  n) 
int  n; 

Passgraph  graphl; 

/•  this  procedure  prints  the  graph  (incidence  matrix)  for  graphl 
which  will  be  a  nxn  matrix  of  0's  and  l's  •/ 

\ 

Int  i.j; 

char  ‘space; 


printf("\n"); 

printf(”  "); 

if  (  (n  +  4)  >  (SCREENSZ/2)  ) 


space  = 


else 

i 


l 


space  =  "  "; 
printf(” 


1 


for  (i=  10;  i  <=  n;  i++)  *  - 

printf("%s%ld'\  space,  i/ 10); 
printf("\n"); 
for  (i=l;  i  <=  n;  i++) 
printf(”%s%ld'\  space,  i  MOD  10); 
printfCNn"): 
for  (i=0;  i  <  n;  i++) 

printf("\n"); 
for  0=0;  j  <  n;  j++) 

printf("7Ss%ld'',  space,  (graphl[i][j/(intsrze)]  &  place[j3)  ?  1  :  0); 
printf("  7.2 d",  i+1); 

printf("\n''); 


. . * . . . 

startgraph(graphl,  n) 

Passgraph  graphl; 

int  *n; 

/•  this  procedure  generates  a  starting  graph  with  *n  nodes. 
The  graph  is  a  ring  with  opposite  sides  attached.  */ 

[ 

int  i.  j; 

Int  distance; 
int  nextpt; 

for  (  ; ;  ) 

prlntfC'Please  enter  the  desired  degree  of  the  graph:  "); 


search.c 

...print  control 


print  graph 


startgraph 
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... startgraph 

3canf("%d'',  ^degree); 
printf("\n"); 

if  ((degree  <  2)  ||  (degree  >  Maxsize  —  1)) 

i 

printf(”The  degree  must  be  an  integer  between  2  and  %d\n",  Maxsize  -  l); 

continue; 

1 

printf("please  enter  the  number  of  nodes  you  want:"); 

scanfO'/Sd”.  n); 

printf("\n"); 

if  ((*n  <=  degree)  ||  (*n  >  Maxsize)) 

( 

printf("The  number  of  nodes  must  be  between  % d  and  %d\n", 
degree  +  1,  Maxsize); 

continue; 

I 

if  ((degree%2  ==  1)  &&  (*n%2  ==  1)) 

printf("Far  odd  degree  the  number  of  nodes  must  be  evenAn"); 

continue; 

1 

else 

break; 

I 

for  (i=0;  i<*n;  i++) 

for  (j=0;  i<=((*n)-l)/(Intsize);  j++)  ‘ 

graph  1  [i][j]  =  OL;  1 


/•  this  section  sets  up  the  ring  •  / 
graphl[0][0]  =  place[0]  ; 

graphl[oi[((*n)-l)/(Intsize)l  |=  place[((*n)-l)]; 
graphl[((*n)-l)][0]  |=  place[0]; 
for  (i=l;  i<*n;  i++) 

graphl 


graph 1 
graphl 


ilfi/(Intsize)]  |=  place[i]; 
iJL(i-l)/0ntsize)]  |=  place[(i-l)]; 
i-l][i/(fntsize)]  |=  place[i]; 


/*  this  section  makes  the  connections  across  the  ring  •/ 
for  (i=l;  i  <=  (degree  —  l)/2;  i++) 

l 

distance  =  1  +  ((i  •  ((*n)  —2)) /(degree  -  1)); 
for  0=0;  j  <  *n;  j++) 

nextpt  =  (j  +  distance)%*n; 
graphl[j][nextpt/(Intsize)]  |=  place[nextpt]; 
graphl[nextpt][j/(Int3ize)]  |=  placefj]; 


l 
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^include  "declarlh" 

diamin£c(vectl ,  diaml)  d.'iCLTTL'i'n.f G 

int  ‘vedtl,  diaml; 

i 

ini  i; 

printf("%ld  improvements  so  far.  ",  status. graphcount); 
printf'("Diameter  vector  for  current  min:\n"); 
for(i=l;  i<=diaml;  i-^+) 
priritf("%5d",ii); 
printf.("\n?); 
for(i=l;  i<=diaml;  i*+) 
pridtf("%5d'V‘(vectl  ■+■  i)); 
print^("\nT); 


takepowers(graphl,  vectorl,  n,  diameter,  olddiam)  tdhspCVJBTS 

Passgrajph  graphil; 

int  vector ltO; 

int  n; 

int  ‘diameter; 
int  olddiam; 

/•  this  procedure' finds  all  'powers'  of  graphl,  which  has  n  pts, 
and  formst  vector  l  which  has  the  count  of  the  number  of  0's 
■in  the  corresponding  power  of  graph  1. 

The  procedure  stops  irJring  powers 

when  two  entries  for  vector!  are  the  same,  or  when  the  number  of 
0's  is  ‘,0,  or  when  olddiam  is  seen  to  bo  at  least  as  small  as  the  new 
diameter,  and  returns  that  as  diameter. 

•/ 

l 

intp: 

int  numzeroes; 
int  cvencdd; 

eveno'dd  =1 0; 

for{p=l;  Cp<alddiam)&&(vectarl[p]  !=  vectorl[p— l])&&(vectorl[p]);p++) 

t 

evenodd==  1  —  evenodd; 

if  (p  ==  a) 

(*multjply)(grapkl.  graphl,  graph[evenodd],  a,  &numzeroes); 

else 

•(♦muRtplyK'graphl ,  graph{  1  -evenc-dd],  grnph{-evenodd],  n,  Anumzeroes); 
vecfcorl[p  IQ  =  numzeroes; 

I 

‘diameter—  p; 


. . . . ... . . . . . . ♦ . . 

takepo(graph! .  vectorl,  n,  diameter,  oidaiam)  takepo 

Possgroph  graphl; 
int  vecborlf}, 
int  n, 

int  ‘diameter, 
int  olddiam;  , 

/•  this  procedure  finds  all  'powers’  of  graph!,  which  has  n  pts, 
and  forms, vector  ]  which  has  the  count  of  the  number  of  0's 
in  thci corresponding  power  of  graphl. 

The  procedure  stops  taking  pov.'crc 
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when  two  entries  for  vector  1  are  the  same,  or  when  the  number  of 
0's  is  0,  or  whan  olddiam  is  seen  to  be  at  least  as  small  as  the  new 
diameter,  and  returns  that  as  diameter. 

•/ 

l 

ini  p; 

int  numzerces; 
ini  evenodd; 

evenodd  =  0, 

for  (p=  1 ;  (p  <  olddiam)  &&  (vectorl[p]  !=  vectorl[p  -l])<3£&(vectcrl[p]),  p-- 

f 

everodd  =  1  — everodd; 
if  (p  ==  1) 

(*multiply)(graph  1,  graphl,  graph[evenodd];  r.,  icnumzerces); 

else 

(*multiply)(graphl,  graph[l— evenodd],  graph[evencdd],  r.,  &r.unizorces); 
printf("\r."}; 

print.gr  a  ph(graph[evenodd],.numpts); 
vector l[p  -  l]  =  numzeroes; 

J 

♦diameter  =  p; 

i 


/*****«.*«********«*****•***••*»**«.******»»»*»»««»•»*«  t  •**»*«**»»««»*««•*/ 

zeroccur.t(graphl(  n,  count) 

Passgraph  graphl; 

int  n; 

int  ♦count; 

/♦  this  counts  the  number  of  0's  in  the  matrix  for  graphl,  and  returns 
that  number  in  count.  Ii  assumes  the  matrix  is  symmetric  and  has  all 
l’s  on  the  main  diagonal.  ♦/ 

l 

int  i.j; 

•count  =  0; 
for  (i=C;  i<r.;  i-r-) 
for  (j=i-l;  j<r.;  j~~) 
if  (  !(graphl[i][j/(Intsize)]'&  placefj])) 

(♦count)-^--t-; 

(♦count)  *=  2; 

i 

/*♦***♦♦ . *♦**♦***♦♦ . . . ♦♦**♦**♦***♦♦♦♦ . ♦•♦•♦ . . 

multgraph(graph  1,  graph.2,;graph3,  n,  zerocr.t) 

Passgraph  graphl,  graph.2,  graph3; 

int  r.; 

int  ♦zerocr.t; 

/*  this  procedure  'multiplies'  the  matrices  for  graphl  and  graph2,  and 
returns  the  product  in  graph3.  All  matrices,  are  symmetric  n  x  n.  */ 

I 

int  i,  j,  k; 

int  numfields  ; 

long  *grpb3,  ♦grphl; 

Jong  placei; 
int  ioverlntsize; 

(•zeroont)  =  0; 
numfields  =  (n-l)/(Intsize); 
for  (i=0;  i<n,  i+-t-) 

grphl  =  ’(graph  1  +  j); 


...takepo 


) 


zerocount 


mult  graph 
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..  multgraph 


grph.3  =  «(graph3  -r  i); 

placci  =  placefi]; 

ioverlntsize  =  i/(Intsize); 

for  (j=0;  j  <=  numfields;  j++) 

♦(grph3  +  j)  =  0L; 

♦(grph3  +  ioverlntsize)  =  placei; 

for  (j=0;  j  <  i;  j++) 

< 

% 

for  (k=0;  k  <=  numfields;  k++) 

(  *(grphl  +  k)  &  grapii2[j][kj) 

1 

*(grph3  +  j/(Intsize))  |=  place[j]; 
graph.3[j][ioverIntsize]  |=  placei; 
goto  nextj; 

(♦zerocnt) +-r; 

nextj: 


i(*zerocnt)  ♦=  2; 


multlgrapk(graphl,  graphs,  graphs,  n,  zerocnt) 

Fassgraph  graphl,  graphs,  graph3; 

iiit  n; 

int  *zerocnt; 

/*.  inis  is  the  same  as  'multgraph'  except  that  it's  'tuned  for  n<=  (ihtsizc)  ♦/ 

i 

lint  i.  j; 

Hong  *grph3,  grp  hi; 

•long  placei; 


mult  1  graph 


l(*zerccnt)  =  0; 

•*(*(graph3))  =  JL; 

Ifor  (1=1;  i<n;  I++) 

w 

grphl  =  *(*(graphl  +  i)); 
grphS  =  *(graph3  +  j); 
(•grph3)  =  placei  =  place[ij; 
for  (j=0;  j  <  i;  j^+O 
if  (grphl  2c  *(graph2[j])) 
f 

*(grph3)  1=  place[j]; 
*(graph3[j])  1=  placei; 

I 

else 

(*zerccnt)++; 

;(*zcrocnt)  *=  2; 


muU2graph(graphl,  graph2,  graph3,  n,  zerocnt) 

Pas3graph  graphl,  graph2,  graph3; 
int  n; 

int  *zerccnt; 

/•'this  is  the  same  as  ' multgraph '  except  that  it's  tvr.cd  for  n  >  (fntsiza)  •/ 


muitSgraph 
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...nrvj.li2gra/ph 

/*  and  n  <=  2  *  Intsize  */ 

l 

int  i,  j; 

long  *grph3,  ‘grphl; 
long  piacei; 

(♦zerocnt)  =  0; 
grph3  =  ‘graphs; 

(*grph3)  =  1L; 

*(grph3  1)  =  0L; 
for  (1=1;  idntsize;  i++) 

1 

grp  hi  =  ‘(graph  1  +  i); 
grph3  =  *(graph3  ■»  i); 

*(grph3  4  1)  =  0L; 

*(grph3)  =  piacei  =  place[i]; 
for(j=C;  j  <  i;  j4+) 

if  ((‘grphl  &  *(gxaph2[j]))  |]  (‘(grphl  +  1)  if.  gr?.ph:2[j][l])) 

*(grph3)  !=  place[j]; 

*(graph3[j])  j=  piacei; 

I 

else 

(•zerocnt)++; 

\ 

for  (i=fr:tsi'je;  i  <  n;  i+4-) 

f 

( 

grphl  =  ‘(graphl  +  i); 
grpr.3  =  ‘(graph.3  4  i); 

*(grph3;  =  0L; 

•(gvph3  -  1)  =  piseej  -  placc-fij; 
for  (j=C;  j  <  i.  j--) 

if  ((*gi'phl  i:  *(.graph2[ij))  |j  (‘(grphl  tl)4  graphi:[j][ !})) 

/ 

I 

‘(grphS  -*•  j/(I.ntsise))  !=  place[j]; 

^  graphS [j](l]  j-  piacei; 

else 

(•zerocr.t)  4+; 

l 

(‘zerocnt;  *=  2; 

J 
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